Density-Functional-Based Determination of Vibrational Polarizabilities in Molecules 
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We develop a direct derivation for the primary contribution to the vibrational polarizability for 
molecules, clusters and other finite systems. The vibrational polarizability is then calculated within 
the generalized gradient approximation to the density functional theory for a variety of molecules 
and clusters. The agreement between theory and experiment is quite good. The results show that for 
small ionic molecules and clusters, inclusion of the vibrational polarizability is necessary to achieve 
quantitative accuracy. 
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I. INTRODUCTION 

Materials with high dielectric constants have many 
important technological applications. For molecular as- 
sembled materials, the Clausius-Mossotti relation tells us 
that the dielectric constant is related in a simple way to 
the molecular polarizability and points to critical densi- 
ties at which this relation is expected to diverge. Predict- 
ing the total second-order molecular polarizability is one 
capability that is required for the computational design 
of materials with high dielectric constants. The second 
largest contribution to the static second-order polariz- 
ability tensor is generally accepted to be due to field- 
induced atomic relaxation. As the effect is of interest to 
several fields of research, a common terminology is lack- 
ing. This effect has been referred to as the displacement-, 
atomic-, nuclear-, relaxation- or vibrational polarizabil- 
ity. Here we adopt the latter term and determine this 
effect within the double harmonic approximation— As 
discussed below, this effect is governed by the dynamic 
effective charge tensor which is known to account for in- 
frared intensities of vibrational modes in molecules and 
clusters. 

The dynamic effective charge tensor describes how the 
total dipole moment of a molecule or other finite system 
changes due to an atomic displacement. For a simple 
dipole consisting of two point charges ±Q, the change of 
the dipole per unit change of the separation is just Q. 
For crystals, the local dipole is not necessarily a well- 
defined notion, and effective charge tensors come in mul- 
tiple forms 2 . However, the lowest nonzero moment of a 
finite system is a well defined quantity. For a neutral 
molecule with N atoms, both the electrical dipole mo- 
ment p and the derivative of the dipole moment with 
respect to the i'th atomic position Ui are well defined. 
Only the latter quantity is uniquely defined for a charged 
molecule. These derivatives may be expressed as a 3 x 3N 



tensor Z which has units of charge, and is written 



Z a ,ifi = dpa/du 



2/i • 



(1) 



This tensor is also sometimes called the "polar tensor"—, 
and it is used for the calculation of the infrared 
intensity&iiii. The infrared intensities are also related 
to the vibrational component of the dc molecular polar- 
izability^, and a simple proof of this is included below. 

By vibrational polarizability, we refer to the follow- 
ing physics. When a molecule is placed in a static elec- 
tric field, it can lower its energy through several mecha- 
nisms. First, the electronic clouds rearrange themselves 
in response to the field which leads to an induced elec- 
tronic dipole moment given by p e \ x = ^yOt c i t xyEy. This 
is generally the largest linear effect. Second, this in- 
duced dipole moment is further modified since the atomic 
positions rearrange themselves in response to the forces 
caused by the direct application of a field and the sub- 
sequent electronic rearrangement. The tensor that de- 
scribes the portion of the induced dipole moment due to 
atomic rearrangement is what we refer to as the vibra- 
tional polarizability. 

To concentrate on effects due to vibrational polarizabil- 
ity, we ignore molecular rotation and assume the molecule 
to be oriented in the lab frame. Equivalently, we work 
in a frame which is tied to the molecule, so that there 
is a fixed dipole moment. The polarizability a is then a 
3x3 tensor which reduces to a scalar for symmetrical 
molecules such as CH4 or SF@. The molecular vibrations 
within the harmonic approximation correspond to the 
classical normal modes of a coupled system of oscillators 



(2) 



where the 3N x ZN force constant tensor K is defined as 
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K, 



d 2 £ 



(3) 



and £ is the total energy of the molecule at zero field. 

The dipole moment p is a first derivative of the energy 
{£ ) and the dynamical charge tensor and the electronic 
polarizability tensor are second partial derivatives given 
by: 

Pa = -d£/dE a , (4) 



Z aM = -d 2 £/dE a du ifi = dF ifi /dE a , (5) 



aw.a/9 = -d 2 £/dE a dE p . (6) 

In Eqs. (4-6), the electronic degrees of freedom must be 
relaxed in response to changes of the independent vari- 
ables (E, u\, un) and the derivatives are evaluated at 
zero field and displacement (E = Ui = 0). Eq. [5] also 
shows that the dynamical charge tensor determines how 
the Helhnann-Feynman (HF) force (F^ = —d£/dui^) 
changes due to the application of an external electric 
field. As discussed in Ref. 0, the relationship between 
the dynamic effective charge tensor and the derivative of 
the HF force is both instructive and optimal for efficient 
determination of infrared and Raman intensities. 

Now the total energy of the molecule may be expanded 
as a Taylor series in powers of both the atomic displace- 
ments and applied electric fields according to 

£ = £ -p ■ E - ^E ■ a e i ■ E - E ■ Z • u+ • K • u. (7) 

In the above equation p and a are the zero-field values of 
the dipole moment and polarizability respectively. The 
tensor notation is fairly obvious except perhaps for the 
asymmetrical tensor Z whose transpose Z T is defined by 



i? • Z • u = u • Z ■ E. 



(8) 



Now if a static external field E is applied, the atomic 
coordinates u will relax to new positions to minimize the 
energy according to 



u = K _1 • Z T • E 



(9) 



and the corresponding energy of the relaxed molecule is 
£ = £ -p-E--E-[a e i + a vih ].-E (10) 

In the above, the vibrational part of the polarizability is 
given by 



"vib = Z • K 



-l yr 



(11) 



In fully indexed Cartesian form, the polarizability matrix 
is 



(12) 



While the above expression clearly exhibits the isotopic 
independence of this part of the polarizability tensor, a 
simpler expression, directly comparable to experimen- 
tal observables, is possible by rewriting this energy in 
terms of the normal modes of vibration. Let \v > denote 
the eigenvector and u> v the corresponding eigenfrcquency, 
which satisfies the Newtonian equations 



K\v >= luIM\v >, 



(13) 



where M is the mass tensor which in the atom displace- 
ment basis (in) is 

Mip ljv = MiSijS^. (14) 
The orthogonality and completeness relations are 

< v\M\v' >= <W (15) 



El- 



>< V = 



(16) 



The force constant matrix can be written as 

K = 5^M|« > uo 2 v < v\M (17) 



J2\v>^ v 2 <v\ 



(18) 



The effective charge tensor can now be written in the 
eigenvector basis as the charge vector for each normal 
mode 



(3-Z\v> 



(19) 



Then the vibrational polarizability can be written as a 
sum of contributions from the normal modes, 



OVib, a/3 — 2. Za,vU 



-2 7 T 



(20) 



This equation is a generalization of a known relation 6 
between infrared intensities and static polarizability. In 
the past, this equation has been used to determine 
vibrational polarizabilities from experimental IR data 
and from calculations^ We include our derivation here 
because it appears to be rather simple in compari- 
son to previous derivations that appear in the litera- 
ture. Eq. (20) follows immediately from Eq. (17.29) of 
Born and Huangpii It has also been derived by Flytza- 
nisii^ Probably the earliest modern discussion of vibra- 
tional polarizabilities using quantum-mechanical deriva- 
tions can be found in Ref. lj]|8j where applications to 
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CHCI3 and CHF3 are discussed and the above formula 
is derived within a sum over states method within the 
clamped nucleus approximation. Eq. (1) and A5 of 
Ref.0b lead to our Eq. (20). However, as noted in Refs 7 
and0, one of the earliest discussions dates back to 1924. 9 
In addition to the interaction discussed above, there 
are other smaller vibrational effects that modify the po- 
larizability of a molecule. The presence of the field mod- 
ifies the spring constant matrix which changes the zero- 
point energy of the molecule. Also, the occurrence of 
anharmonicity, both diagonal and off-diagonal, leads to 
further corrections. We are unaware of discussions on 
the role of off-diagonal anharmonicity, but discussion of 
the zero-point effect and diagonal anharmonicity may be 
found in Ref. I n the notation of the work of Marti 
and Bishop, the above term is equivalent to [/i 2 ] 0,0 in 
their paper. 



II. COMPUTATIONAL DETAILS 

The calculations presented below have been performed 
using the NRLMOL suite of density-functional-based 
cluster codes^ The Perdew-Burkc-Ernzcrhof (PBE) en- 
ergy functional has been used in all calculations^ The 
Kohn-Sham equations are solved self-consistently for 
each electron in the problem. Then the HF forces are 
calculated and the geometries are updated using stan- 
dard force optimization methods. Geometries were con- 
sidered converged when the force on each atom fell below 
0.001 Hartree/Bohr. However, for the Na and H2O clus- 
ters we used a tighter force convergence criteria of 0.0001 
Hartree/Bohr. The numerical integration mesh was also 
significantly more dense for our calculations on the wa- 
ter molecules. The method for generating the basis sets 
used for these calculations is discussed in Ref. EE Basis 
sets and the unpublished geometries are available upon 
request. Once the optimized geometriesiii^iiS are ob- 
tained, the vibrational frequencies, eigenvectors, and dy- 
namical charge tensors (Z™^) are determined using the 
method discussed in Ref. 5f. We then use Eqn. (20) to de- 
termine the vibrational component of the polarizability. 
As discussed in Ref. 0, the infrared and Raman spectra 
showed some sensitivity to the inclusion of longer range 
polarization functions. We have used such polarization 
functions for the calculations displayed in Table. I. 



III. RESULTS 

Table I presents calculations on several molecular sys- 
tems which include both covalent, ionic bonding and hy- 
drogen bonding. It also include calculations on systems 
with both loosely and tightly bound electrons. We have 
calculated both the electronic and vibrational contribu- 
tions to the polarizability tensor. Agreement is generally 
good. 



Fullerene Molecule: The polarizability of the 
fullerene molecules has been well studied both theoret- 
ically and experimentally. 20 ' 21,22 23 ' 24 Here we calculate 
the electronic polarizability to be 82.9 A 3 which is in good 
agreement with one of the earliest density-functional cal- 
culations^ - of 83.5 A 3 . This earlier calculation used the 
same code, a slightly different version of DFT, slightly 
smaller basis sets and geometries that were not as well 
converged. 20 The good agreement between the early and 
most recent calculations indicate that the electronic part 
of the neutral fullerene polarizability is rather robust, 
and the experimental polarizabilitySisSSiSi 2 ^ is known to 
be very close to this number as shown in Table I. Based 
on experiments, it has also been suggested that the po- 
larizability due to lattice relaxation is 2 A 3 — which is 
small but still four times larger than the value calculated 
here. The deviation may be due to the lower T^ sym- 
metry that occurs when the icosahedral Ceo molecules 
are placed on a cubic lattice. Such a symmetry lower- 
ing would cause some of the optically silent G„ and H u 
modes to split and partially fall into the IR active Ti u 
manifold which in turn could lead to additional vibra- 
tional polarizability. There will also be weak IR activ- 
ity due to weak intermolecular vibrations activated by 
weakly broken translational symmetry^ Also included 
in Table I are the electronic and vibrational polarizabili- 
ties of a Ceo molecule with an endohcdral Kr atom. The 
addition of the Kr atom adds another infrared mode due 
to a rattling motion of Kr inside the Ceo cage. The low 
frequency Kr rattling mode is found to be at 88 cm -1 
but the IR intensity associated with this mode is 1000 
times smaller than the four Ti„ modes associated with 
the fullerene cage. Because of this the vibrational polar- 
izabilities are unchanged due to the addition of an inert 
endohcdral atom. 

Acetylene: The acetylene molecule provides an inter- 
esting test case because the anisotropy of the polarizabil- 
ity tensor is reversed significantly by the inclusion of the 
vibrational terms. For example, in A 3 the electronic and 
vibrational polarizability tensors have been measured to 
be (2.43, 2.43, 5.12) and (0.667, 0.667, 0.027), respec- 
tively. Density functional theory yields (2.96, 2.96, 4.78) 
and (0.71, 0.71, 0.030) A 3 which is in reasonably good 
agreement with experiment. 

Halogen containing Ionic Molecules: Halogen 
containing compounds are known to exhibit high vibra- 
tional polarizabilities as would be expected since they 
make very good ionic systemsAS We have performed cal- 
culations on NaF, SiF6, SiF4 and TiC^. Of the molecules 
in this size regime listed in the large database of Gussoni, 
the latter three stand out as having very large vibrational 
contributions. The agreement between theory and ex- 
periment is in the neighborhood of 15 percent for these 
systems. 

Isomeric dependencies: Acctonitrile (CH3CN) and 
methylisonitrile (CH3NC) have the same chemical com- 
position. However, the former has the two carbon atoms 
bound to one another while the latter has the nitrogen 
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bound to the methyl radical. This causes a five per- 
cent difference in the electronic polarizability and a fac- 
tor of two difference in the vibrational polarizability. The 
source of the deviation in the vibrational polarizability is 
clearly due to changes in the spring constant matrix since 
Eq. (11) shows that changes of mass cannot perturb the 
vibrational contributions within the double harmonic ap- 
proximation. 1 There is a one-fold mode at 2269 1/cm for 
acetonitrile that is reduced to 2149 1/cm for methylisoni- 
trile. In addition to a reduction in the vibrational fre- 
quency, the infrared intensity of the methylisonitrile is 
2.57 compared to 0.227 in the case of acetonitrile. This 
mode accounts for about 75 percent of the difference in 
the vibrational polarizability. The large change in in- 
frared intensity in this frequency range should be a clear 
indicator of methylisonitrile isomerization to acetonitrile 
at higher temperatures. 

Sodium clusters: In two recent papers 2 ^ Blundcll, 
Guet and Zope and Kronik, Vasiliev and Chelikowsky 
have calculated the temperature dependence of polariz- 
abilities in sodium clusters. They show that temperature 
effects enhance the apparent polarizability at 300K. This 
temperature enhancement appears to account for most 
of the difference between experiment and the calculated 
electronic polarizabilities from many different theoreti- 
cal calculationsi 29 i 30 i 31 i 32 i 33 i 34 Our results show that the 
vibrational contribution to the polarizability is indeed 
small for the sodium clusters which supports the asser- 
tion that temperature effects are important in these sys- 
tems. 

Weakly Bound Molecules: As mentioned in our 
discussion of fullerene molecules the vibrational polar- 
izability between two weakly bound molecules could be 
enhanced if the weak intermolecular vibrations are IR ac- 
tive. As discussed in Ref. |35| the water dimer represents 
an extreme example of this case. As shown in Table 
I, we find the electronic polarizability of this molecule 
(3.19 A 3 ) to be approximately twice that of a water 
monomer. The electronic polarizabilities obtained for the 
water trimer and pentamer also show a linear scaling as 
a function of the number of molecules. This result is 
in good aggreement with the work of Maroulis et 
(2.90 A 3 ) and Eckart et al (2.48 A 3 ), 35 Maroulis et 
have carefully studied the electronic polarizability as a 
function of both basis set and level of correlation. The 
uncertainties due to these effects are at most 12.5 per 
cent indicating that large deviations from these values 
must be due to other effects. Our calculated double- 
harmonic vibrational polarizability of 23.30 A 3 is indeed 



a factor of seven times larger than the electronic polar- 
izability. Eckart et al find this term to be even larger 
(39.2 A 3 ) and further demonstrate that anharmonic cor- 
rections enhance the vibrational component of the dimer 
by an additional factor of three. The large vibrational en- 
hancement in polarizability in going from the monomer 
to dimer is indeed interesting. In particular, the scal- 
ing of this term as a function of system size is impos- 
sible to guess based upon the results of the monomer 
and dimer. It is reasonable to expect that this large re- 
sult should be an upper limit since a dielectric medium 
that is coupled to an IR active mode should counteract 
the IR activity and thus the vibrational polarizability. 
So the presence of more water molecules should lead to 
a vibrational polarizability that is eventually sublinear 
in the total number of molecules. To partially address 
this point we have performed additional calculations on 
the trimer and pentamer. Our results show a decrease 
in the total vibrational polarizability in going from the 
dimer to the trimer and a flattening of the total vibra- 
tional polarizability for the pentamer. For the pentamer 
the ratio of the vibrational to electronic polarizability 
has decreased significantly from seven for the dimer to 
slightly less than two for the pentamer. Overall, these 
results show that weak intermolecular vibrations can en- 
hance the vibrational polarizability over what is deter- 
mined from intramolecular vibrations. 



IV. SUMMARY 

We have presented a straightforward derivation for the 
vibrational polarizability of a molecule. We have used the 
generalized-gradient approximation to the density func- 
tional theory to evaluate both the electronic polarizabil- 
ity and this vibrational correction. In accord with exper- 
iment, our results show that this term can be important 
in smaller ionic molecules and weakly bound systems, but 
that it is smaller in covalent systems or where the frontier 
electrons are delocalized. 
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TABLE I: Calculated and experimental vibrational polariz- 
abilities (A 3 ) for molecules. This is one third of the trace of 
the polarizability tensor. Unless otherwise stated the exper- 
imental values are taken from Ref. 0- Experimental data for 
Cgo is from Refs. |2lll22l|23!l24| and references therein. For the 
pure sodium clusters, the experimental values are total po- 
larizabilitics (denoted with a *) which have been taken from 
Refs.|26j27l 



Molecule 


Vibrational 


Electronic 




Theory 


Expt. 


Theory 


Expt. 


H 2 


0.04 


0.04 


1.57 


1.45 


(H 2 0) 2 


23.30 




3.19 




(H 2 0) 3 


13.50 




4.82 




(H 2 0) 5 


14.52 




8.13 




NF 3 


1.15 


0.70 


3.07 


2.81 


HCCH 


0.48 


0.45 


3.56 


3.40 


CH 4 


0.04 


0.03 


2.63 


2.60 


TiCk 


2.04 


>1.68 


15.03 


15.0 


SF 6 


2.29 


2.29 


5.15 


4.49 


SiF 4 


2.09 


1.75 


3.72 


3.32 


HCN 


0.21 


0.14 


2.62 


2.59 


H 3 CNC 


0.07 




4.87 




H 3 CCN 


0.03 


0.04 


4.59 


4.28 


C 60 @Kr a 


0.55 




83.3 




C60 


0.58 


2.0 


82.9 


83 


Na 2 


0.00 




35.91 


37.91[*] 


Na 3 


1.72 




60.89 


69.8[*] 


Na 8 


0.63 




116.2 


133.5[*] 


NaF 


0.29 




2.71 




TigC5 2 


3.72 




56.40 




As@Nii 2 @As 20 


4.56 




140.86 





c Ref. Il7ll8ll9l respectively 



